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Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

A hydrodynamic model is postulated for a small, liquid or solid sphere 
film boiling on a saturated cryogenic liquid surface, such as may occur in 
propellant spillage accidents and in the preparation and preservation of 
some biological species (e.g. , blood and enzymes), The momentum and 
energy equations are solved analytically to predict the velocity and temper- 
ature fields in the vapor film beneath the sphere „ This leads to a heat- 
transfer coefficient of the form 
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k 

R 




2 k Pd pgX * 
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where k is the thermal conductivity of the vapor, p d the density of the 
sphere, p the density of the vapor, X* the modified latent heat, R the 
radius of the sphere, p. the vapor viscosity, T^ and T g the temperature 
of the sphere and the liquid, respectively, and f{0*) a function which 
depends on the characteristics of the sphere and the supporting fluido The 
values f(6*) and dimensionless forms of the heat -transfer solutions are 
given in the paper. This expression for the heat-transfer coefficient com- 
bined with a thermal balance predicts the time a hot sphere will remain 



levitated in Leidenfrost film boiling when placed on a cryogenic liquid. The 
agreement between experimental and predicted times is fair, 

A motion picture supplement demonstrating the phenomenon is available 
upon request from the authors, 

INTRODUCTION 

Film-boiling heat transfer between two liquids at different temperatures 
has been under extensive investigation recently in consideration of propellant 
spillage accidents (ref, 1), If liquid propellants spill accidently during test- 
stand or launching operations, a catastrophic expolsion resulting from the 
detonation of the fuel and oxidizer is possible. In particular , a number of 
experiments (ref, 1) have been performed involving the mixing of small kero- 
sene spheres in liquid oxygen. When small kerosene spheres are dispersed 
in a matrix of liquid oxygen, the sensible heat of the kerosene spheres is 
drained off and vaporizes the liquid oxygen, thereby forming a film of oxygen 
vapor around the sphere. In film boiling, the vaporization rates are one or 
more orders of magnitude less than those of peak-nucleate boiling. In nucleate 
boiling, extreme turbulence and mixing will occur because of the large amount 
of vapor sporadically generated at the interface. On the other hand, film 
boiling is a relatively quiescent phenomena, and as such can have a significant 
damping effect on the mixtures explosive potential. The encapsulation of the 
liquid sphere (as it freezes), as shown herein, can also provide significant 
dampening. 

The liquid-liquid film boiling phenomenon is also of interest purely from 
stability considerations. In this case, a liquid sphere with a higher specific 
gravity floats upon a cryogenic liquid of lower specific gravity while film 
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boiling is taking place . The sphere usually freezes and continues to float 
(solid- liquid film boiling) then suddenly at the transition from film boiling 
to nucleate boiling, the vapor blanket surrounding the sphere disappears. 
Then the sphere sinks beneath the liquid and falls to the bottom of the pool. 

This paper considers two aspects of this liquid-liquid or solid-liquid 
film-boiling problem. First, a hydrodynamic model is postulated for a 
small, spherical, liquid drop film boiling on a cryogenic surface. The 
conservation equations of momentum and energy are solved analytically 
to predict the velocity and temperature fields, and a simple theoretical 
correlation for the heat-transfer coefficient is presented. Details of the 
solution are presented in reference 4. 

Secondly, a heat balance is performed on the sphere, which leads to 
a correlation equation for predicting the time required for a liquid sphere 
to cool from room temperature to its freezing point and then to the temper- 
ature at which boiling ceases; liquid subcooling was neglected. Theory is 
checked with experiment by placing water spheres on liquid nitrogen and 
comparing the experimental and theoretically predicted times for the 
sphere to "float”. Fair agreement was found. 

A motion picture supplement demonstrating the phenomena is available 
from the authors upon request. 

BASIC MODEL AND EQUATIONS 

Consider the film model depicted in figure 1. This model appears to 
fit the physical case of a solid or liquid sphere ’’floating” on a second liquid 
(solid- liquid and liquid- liquid film boiling). For example, drops of water 
floating on a sea of liquid nitrogen. There is, of course, no reason to limit 



4 


the model to liquid spheres or cryogenic fluids. Slush and small solid 
materials could also be utilized on any liquid provided the proper temper- 
ature difference and buoyancy criteria are maintained. 

The model may be applied either to the case where the sensible heat 
of the floating sphere evaporates the encompassing liquid or where the 
sensible heat of the encompassing liquid evaporates the floating sphere, 
(e.g. , liquid nitrogen sphere on water). For simplicity, it is assumed 
throughout the remainder of this paper that the supporting fluid evaporates 
at the interface. Physically, for the experimental results considered 
herein, the heat passes from the warmer spherical drop and evaporates the 
fluid material beneath it, thereby forming a supporting vapor layer under 
the sphere. The vapor layer is, in turn, supported by the interface. Sur- 
face tension, density ratio, buoyancy, surface curvature, and the ” non- 
wetting’ ' character of the interface play an important role in determining 
whether the interface is strong enough to support the sphere. During the 
evaporative, or film-boiling, process the temperature of the sphere de- 
creases until it coincides with the Leidenfrost temperature of the liquid. 

At this time film boiling ceases, and transition boiling quickly followed by 
nucleate boiling begins. The vapor layer beneath the sphere no longer 
exists; rather, small nucleate bubbles form in the sites on the solid (liquid) 
sphere. Thus, because there is no longer a supporting vapor gap, the 
sphere sinks beneath the surface and falls to the bottom. 

The following assumptions are made in developing the model: 

(1) Evaporation and internal circulation of the liquid sphere and radia- 
tion are considered small. 
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(2) The model has complete symmetry with respect to the <& coordin- 
ate. 

(3) The flow of vapor is considered laminar and incompressible; and 
the inertia terms in the Navier-Stokes equations are neglected. (see 
refs. 2 and 4). 

(4) The velocity and temperature profiles are assumed to be in steady 
state and the vapor gap thickness to be a constant. 

(5) At any instant of time, the sphere is at an average temperature T d , 
and the evaporating liquid is assumed to be at the saturation temperature T . 
The properties of the flow field are evaluated at the film temperature 

Tf = (T d + T )/2 and are considered to be constant. This assumption has 
worked quite well, (e.g., ref. 2). 

(6) The major mode of heat transport is by conduction (see refs. 2 and 4). 
Thus, it is assumed that 

k 

l«^±(v 2 ™\ ( 1 ) 

r 30 r 3r 2 dr \ dr J 

r ' 


(i.e. , the convection terms may be neglected). Furthermore polar conduc- 
tion is less than radial conduction, thus 
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(7) Heat transport within the sphere is by conduction alone. With cool- 
ing at the lower surface there should be no instabilities or ensuing cellular 
motion. Supercooling during freezing is neglected. 
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(8) The supporting liquid is at the saturation temperature. Thus, with 
assumption 7, all the heat reaching the liquid produces vapor. 

The governing equations are as follow (ref. 3): 

1. Momentum 


0 = -f£ + „ (v'V . fli . A . 2 1 ± 


cot 0 


p 3r 


r 2 r 2 30 r 2 


g c 9p (2 2 ^9 
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where 


2. Energy 


V 2 = J. A fr 2 ±\ h- L ± fsin 0 

r 2 3r V 3r / r 2 sin 0 39 V 89 / 


0 = JL fr 23 !' 
3r V 3r , 


3. Continuity 


0 = A A (r 2 V ) + — i — A (Vn sin 0) 
r 2 3r r sin 0 30 


(3) 

(4) 

(5) 

( 6 ) 
(7) 


The boundary conditions on the velocity and temperature fields are 
given in table I, where 0* is defined in figure 1, and 6 is the vapor gap 
thickness. 
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Drop surface 

r = R 
9 = 9 

v r = o 
v e = o 

T = T d 

(8) 

Liquid- vapor 
interface 

r = R + 5 
9 = 9 

v r = V r (R +5,9) 
V 9=° 

T = T S 

(9) 

,, Stagnation M 

region 

R < r < R +6 
0 = 0 

<J 

II 

o 

(10) 

* ’Separation” 
region 

r = R 
0 = 0* 

P = P 0 

(11) 


Table I. - Boundary Conditions for the Floating Sphere, 0 ^ 0* 


The boundary conditions are still incomplete at this point since 6 , 
V r (R +8), and 0* are unknowns. Hence, three additional mathematical 
constraints are necessary to make the problem tractable. 

Static Force Balance 
(Neglecting Forces Internal to Drop) 

One additional constraint requires that the weight of the drop be bal- 
anced by the shear and pressure forces acting beneath the drop. These 
forces are depicted in figure 2. Summing the forces acting in the vertical 
direction gives 


0=w d - 



(P cos 9 + T r0 sin 9 - r rr cos 9 ) 


r=R 


PT sin 9 d0 d$ 



P Q cos 9 R 


sin 9 d0 d<& 


(12) 
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The last term in equation (12) represent ambient pressure force acting on 
the upper portion of the drop. 

Interface Energy Balance 

The second constraint necessary for the solution of equations (3) to (7) 
is the interface energy balance. Because the supporting fluid has been as- 
sumed to be at the saturation temperature (assumption 8), all the heat 
leaving the drop produces evaporation of the fluid material. Mathematically, 
this constraint is expressed as 


-pAV ( R + 6,0) = -k~ 
r 0r 


(13) 


r=R+5 


where A is the latent heat of vaporization and -k(3T/3r)| r _ R4 ^ is the con- 
duction heat flux to the boundary of the supporting liquid. 

Free Surface Pressure Head 

The third and final constraint is that the forces supporting the sphere 
must be transmitted and balanced by the supporting fluid interface. Here 
the free surface head Z Q , as shown in figure 3 must be determined from 
a balance of forces acting on the liquid- vapor surface membrane. The 
determination of Z Q such that all forces are balanced and the sphere is 
supported also uniquely prescribes a value to 9*. 

The solution of the governing equation is long and involved, but quite 
straight forward in that the number of equations now matches the number 
of unknowns. The derivation is given in appendix B. 

The expression for the heat transfer coefficient to the sphere becomes 



R 


2 PP d gk X* 

9jLtR(T d - T g )(cos 2 9 * - 2 cos 9 * + 1) 


1/4 


(14) 


where the value of 9 * is given in figure 4, for the optimum (most massive) 
levitated sphere. A first order estimate of 9 * for non-optimum (lighter 
than optimum) spheres is found in reference 5, or calculated directly from 
the technique in appendix B. 

COMPARISON OF EXPERIMENT AND THEORY 

A simple test of the analytical model is performed using a water sphere 
on liquid nitrogen. At the expense of the internal energy of the water sphe- 
roid, liquid nitrogen is vaporized beneath the drop. This vapor forms the 
supporting cushion for the drop, which appears to "float” on the much colder 
nitrogen surface. 

A slight amount of dye coloring was added to the water. During the 
freezing process, the dye color changes and continues to change until the 
entire drop is frozen. The frozen drop is then a hard sphere "floating" on 
a sea of liquid nitrogen. Its free floating state is sharply terminated at the 
onset of nucleate boiling, and the drop falls beneath the surface. 

The size of the water droplet used can be determined by measuring the 
residue (the hard sphere of ice) or by using a pipette. Both techniques were 
used. A stop clock was used to determine the elapsed time between placing 
the water on the nitrogen and the color change in the drop, the time for the 
drop to be completely frozen, and the time for the drop to fall to the bottom. 
The experimental results are presented in table II. 
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To use the experimentally measured drop freezing times as a check 
on the theoretically determined heat transfer coefficient, it is necessary 
to determine the temperature of the drop as a function of time during 
freezing. There are, however, many facets of the freezing problem that 
we do not understand. For example, the convective heat transfer coeffi- 
cient from the top of the droplet as well as internal convection are unknown 
if supercooling occurs (up to T^ - T ice = 40° K is possible), the solid will 
not form at the freezing temperature which completely alters a heat bal- 
ance; the compressive effects of freezing on the interface freezing temper- 
ature and interface growth rate are not known; the effects of cracks, ice 
"worms", and the 4° C inversion point are also unknown; to what extent 
does eccentric cooling and freezing alter a concentric drop freezing anal- 
ysis?; the effect of a pseudo-steady analysis as compared to a complete 
transient analysis is difficult to assess; to what extent is the metastable 
Leidenfrost phenomenon influencing the results ? If a droplet followed the 
metastable Leidenfrost Line (ref. 6) then floating times would be greatly 
extended. 

Each of these unknowns constitutes a report in itself. As a first order 
approximation to the actual problem, two simplified models will be used. 

The first model, called the pseudo-steady state model, takes into ac- 
count thermal gradients in the freezing sphere. The calculated tempera- 
ture-time history for this model is given by the lower curve shown in fig- 
ure 5. The derivation for this model is given in reference 4. 

The second model, called the Newtonian model, assumes that no 
thermal gradients exist in the drop. This is analogous to the classical 
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Newtonian cooling problem. This model will now be considered. First, 
the time required for the drop to cool to the freezing temperature, t , 
must be calculated. Next, the time required for the heat of fusion to be 
removed, tj , must be calculated. Finally, the time for the drop to 
cool to the transition-nucleate boiling temperature (Leidenfrost point) 
is calculated. 

In calculating t for the Newtonian model, the surface of the 
water droplet is assumed to be at 0° C and to remain at that temperature 
until the remainder of the drop reaches 0° C. The sensible cooling time 
can be found from an energy balance 

* _ m d C p,d< T do- T s> (1 „ 

Sen ■ K.tT.ce - V (15) 

where h is determined by the fraction of the surface in film boiling that 
in natural convection 

h = h^l - cos 9 *) + h.n. c. (1 + cos 9 *) (16) 

Experimentally, the natural convection coefficient, h.n.c. , was determined 
to be 


Nu = h ° n ° c ‘ R = ip (17) 

k 

and the film boiling heat transfer coefficient, h^, is given by equation (14). 
For a 0.6 centimeter diameter water drop in nitrogen, the value of t gen 
is 4.2 seconds, as shown in figure 5. 
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The time required to freeze the entire sphere, isothermally, can be 
calculated in the same manner only C p Q - T g ) is replaced by y, 

the latent heat of fusion, that is 



m * y 

^) d (T ice - T g ) 


(18) 


For the 0.61 centimeter water drop on nitrogen, the value of t. was 

1. Lc 

12. 5 seconds. This is represented by the horizontal line shown in figure 5, 
Once the entire sphere has turned to ice, the solid drop will lose its 
sensible, energy and cool to the temperature of the supporting liquid. The 
sensible cooling time can be found from an energy balance on the drop of 
the form 


h d A(T H - TJ = -(pC),„ a V 


3(T d - T g ) 


ice 


3 T 


(19) 


where 


= t - (t 


sen 


+ ‘ice) 


(20) 


and at t = 0, the surface temperature T(R) is taken as the freezing tem- 
perature T. . The solution of equation (19) for the time t, assuming 

ICG 

that varies little with temperature, becomes 


R(pC). 


,T = 


ice 


In 


3h, 


T d- T s 


T. - T 
ice s 


( 21 ) 


*A good assumption except near the Leidenfrost temperature. 
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Thus, the time for a drop to freeze to some temperature beneath its freezing 
point is 


* ^sen + hce + ' 


R(pC) 


ice 


In 


3h 


T. - T 
ice s 


T d- T s 


( 22 ) 


From this equation, the temperature fall off shown by the upper curve in 
figure 5 is predicted. 

As mentioned earlier, as the free floating drop cools and its outside 
surface temperature approaches the saturation temperature, film boiling 
will cease and transition-nucleate boiling will begin. At the onset of 
transition-nucleate boiling, the drop is wetted and sinks. The tempera- 
ture at which this transition occurs is called the Leidenfrost temperature. 

The Leidenfrost temperature of an ice-nitrogen combination is required. 
At present, there is no available data from which to estimate this Leiden- 
frost temperature. Consequently, an estimate of the ice-nitrogen Leiden- 
frost point is made from the data shown in figure 6. 

First, the dimensionless Leidenfrost temperature difference 0^ for 
smooth or rough surfaces, shown in figure 6, is assumed to vary linearly 
as a function of the natural logarithm of (/3/ /3 gt ) » The j3 represents a mea- 
sure of how easily energy can be transported to a metallic surface. The 
variation of (3 with surface roughness has not yet been assessed and the 
superposition principle with be applied using the surface roughness data of 
Cumo et al. (ref. 7). Cumo et al„ (ref. 7) found the Leidenfrost tempera- 
ture difference to be 27 percent higher on a rough surface as compared to 
a smooth surface. Thus the upper curve of figure 6 represents an estimate 
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of how the effects of surface roughness could alter this Leidenfrost temper- 
ature, The surface of the freezing water drop with its crystalline protru- 
sions would be a rough surface. It must be emphasized that this is only a 
first order estimate. The data and lines shown in figure 6 do not take into 
account property variations of the liquid, only the metallic surface proper- 
ties are considered. For ice, the ratio of j3 ice //3 gt is 13.7. Consequently, 
the Leidenfrost temperature difference for ice, a rough surface, on nitrogen 
will be approximately 2.4 times the value measured for water on steel, or 
about 126 K as is illustrated by the horizontal line labeled T L in figure 5. 
The surface temperature of the 0.61 centimeter water drop falls to this 
value in 29.2 seconds, at which time transition-nucleate boiling begins and 
the drop falls beneath the surface. Thus, the theory is seen to be in rea- 
sonable agreement with the data range illustrated by the cross hatched sec- 
tion of figure 5. The interface configuration used in this calculation is illus- 
trated in figure 7. 

Vaporization Rate 

Continuity at the interface requires that 


pV r A 


^cp 


P ^ 
dt 


Z=cp 


dM 

dt 


(23) 


Changing to dimensionless form and substituting equation (B30) and (B31), 
the vaporization rate may be expressed in terms of the vapor gap 
thickness, A, and 0*; (see also eq. (B44)). 

*Ref. 8 gives a Leidenfrost temperature difference, T-^ - T = 15 K, 
and ref. 9 gives 33 K. Therefore, the corrected Leidenfrost temperature 
could range from 113 to 151 K with times of 34.2 and 23.5 sec (see also 
fig. 5). 
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dm / i f 1\ /1 - cos Q* \ _ ( h I R \/l - cos 9* 

dt \ A/\ 2 / \k )\ 2 


(24) 


where 

m = — — — 

4ttRjuS 

Equation (23) may also be formed from the interface energy balance as 


dM 

dt 


h I A AT 47rR 2 AT /I - cos e*\ 


X 


= (47tRju) 




(25) 


If X is replaced by X*, then equations (25) and (24) are equivalent. For 
the 0,6 centimeter diameter sphere, evaluating properties at 286 K, where 
hj = 0.0155 watt per square centimeter per second, the vaporization rate is 

— = 0,013 gm/sec (26) 

dt 


A crude comparison of the vapor generation rate at peak-nucleate boil- 
ing with film boiling near the Leidenfrost point gives a upper limit to the 
ratio of vaporization rates for our example 



p. n. 


dM \ 
dt ) 


f.b. 


* h AT >p.n, 


(0. 9)(16) 3 _ 37 5 
(0. 016)(40) 2 


(27) 
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and a reasonable lower limit becomes 



0.013 


while the vapor generation rate for film boiling increases with increased 
AT, for our example, it is still about an order of magnitude less than the 
peak-nucleate rate. In addition, the generation- rate in film boiling is quite 
uniform when compared to the sporadic nature of nucleate boiling. 

CONCLUSIONS 

Theoretical expressions for the heat transfer coefficient between a 
liquid drop floating in film boiling on a cryogenic fluid were derived from 
the fundamental equations of momentum and energy. The expression for 
the heat transfer coefficient and the quasi-steady freezing analysis lead to 
a reasonable prediction of the time necessary for a drop at room tempera- 
ture to freeze when placed on the cryogenic fluid. Further analysis of the 
freezing process (which is quite complex) and instrumented spheres is 
required. 

The heat transfer coefficient to a sphere floating on the surface in 
film boiling can be shown to be equal to 



k gp d p^* 

Rp(T d - T g )(cos 2 0* - 2 cos 0* + 1) 


1/4 


The value of 9 * may be approximated using figure 4, for optimum levi 
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tated spheres. Non-optimum sphere estimates of 9 * can be calculated 
using the technique of appendix B» 

Two separate and distinct transitions are seen in the film boiling of a 
room temperature liquid drop on a cryogenic surface. First, the drop 
cools from room temperature and freezes. For water drops containing a 
blue dye, this freezing process becomes visible to the eye by a very 
noticeable change in the dye color. Secondly, when the drop temperature 
cools to the Leidenfrost point the supporting vapor layer is destroyed and 
the drop falls beneath the surface. 

A motion picture supplement is available from the authors upon 
request. 


*The "heaviest 1 * sphere that can be floated on a given fluid. 
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APPENDIX A 

SYMBOLS 

2 

A surface area of drop, cm 

B Q Bond number, [(p^ - p)gR 2 ]/ag c 

Cfce specific heat of ice, joule/(g)(K) 

0^ specific heat at constant pressure, joule/(g)(K) 

C specific heat at constant pressure of drop, joule/ (g)(K) 

Pj u 

f(0*) function, (1-2 cos 0* + cos^ 0 *) 

2 

g local acceleration of gravity, cm/sec 

g gravitational constant in Newton’s law of motion, 1 — — 

u 2 

? (dyne) (sec ) 

h heat transfer coefficient to drop, W/(cni )(K) 

hj heat transfer coefficient to the interface, jouie/(sec)(cm ) 

k thermal conductivity, joule/(sec)(cm)(K) 

L reference length, -J^S c /[lp l - p)g) , cm 

M mass, g 

m dimensionless mass x time, M/4?rRpS, sec 

m^ mass of water drop, g 

2 

P pressure, dyne/cm 

2 

P Q atmospheric pressure, dyne/cm 

R radius of drop, cm 

r radial coordinate, cm 

T temperature, K 

drop temperature, K 
T do initial drop temperature, K 

Tj- film temperature, (T^ + T )/2, K 



T. 


ice 


L,st 


T 

T 


s,H 9 0 


d 


t. 

ice 

t 

sen 

V 


V r 

v e 

V 

W 

X 

Z 

z 

0 
0 . 

^st 

y 

8 

0 L 

0 

0 * 


o 


ice 


freezing temperature, K 
Leidenfrost temperature, K 
Leidenfrost temperature of steel, K 
saturation temperature of liquid, K 

saturation temperature of water, K 

time, sec 

time for heat of fusion to be removed from the drop, sec 
time necessary for sensible heat to be removed from drop, sec 

3 

volume of drop, cm 

radial velocity component, cm/sec 

theta velocity component, radian/sec 

volume, cm ^ 

weight of drop, dyne 

horizontal coordinate, cm 

vertical coordinate, cm 

initial pressure head, cm, fig. 3 

V( k Pm C p)> Qcm 2 )(K)/joule] 2 sec 

/3 evaluated for ice properties, [(cm 2 )(K)/w] sec - '*' 

/3 evaluated for stainless steel properties, [(cm 2 )(K)/w] sec - * 
latent heat of fusion, joule/g 
gap thickness, cm 

< T L - T s>/< T L,st - T s,H 2 0* 
angular coordinate, radians 

defined in fig. 1, radians 
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X* 


P 

v 


P 


Pd 



a 

r r9 ’ T rr 

T 




latent heat of vaporization, joule/g 
modified latent heat of vaporization, joule/g 
dynamic viscosity, g/(cm)(sec) = (poise) 
kinematic viscosity, cm /sec 

O 

density of vapor, g/cm 

O 

density of drop, g/cm 

O 

density of ice, g/cm 

O 

density of liquid, g/cm 

O 

density of metal, g/cm 
surface tension, dyne/ cm 

O 

shear stresses, see fig. 2, and ref. 3, dyne/cm 

time, t = t - (t „ - t. ) 

’ x sen ice ; 

angular coordinate, see fig. 1, radians 
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APPENDIX B 

SOLUTION OF THE GOVERNING EQUATIONS 
Momentum Equations 

The governing equations are made nondimensional to generalize the 
solution by selecting the following parameters 


Re 


c =- 

R 

(Bl) 

V r 

Vir _ 

(B2) 

V 0 

v e = — 
u 

(B3) 

* 

II 

SJ k 

(B4) 

^ = 1 Reynolds number 

(B 5) 



p= p 

(B6) 

*2 

pu* 


e c 


Wd= l T p d R 3 X 
3 s c 

(B7) 

nn * 2 r >2 

w * _ P u R 

(B 8) 

Sc 
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Ra* = 


Pr 


Pr = E 


pip 1 ~ p)gR‘ 

2 C 
/i S 


V T d - T s> 

X* 

Prandtl number 
modified Rayleigh number 


(B9) 


(BIO) 


(Bll) 


Nondimensionalizing the momentum equations (3) and (4) using (Bl) 
through (B6) and introducing the stream function for spherical coordinates 
(ref. 3, p. 131) equations (3) and (4) can be reduced to 


E 4 ^) = 0 


(B12) 


where 


E 2 =i( 


3? 


+ 


sin 6 9 


1 0 


c 


2 36 \sin 6 36 


(B13) 


The similarity transform 


9 ) = F(C)(1 - cos 6) 


(B14) 


is used to convert equation (A 12) into the ordinary differential equation: 


F*(?) = 0 


(B15) 


This is a linear homogeneous equation of fourth order which has a solution 
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F(0 =Cj +C 2 C +C 3 C 2 +C 4 ? 3 


(BX6) 


Therefore, the velocity components may be written as 


V, 


F(0 


v 0 = 


F’(0 tan- 
2 


(B17) 


(B 18 ) 


Applying the boundary conditions (eqs. (8) to (10)), and using equations 
(B16) to (B18), allows the four constants to be expressed in terms of one 
constant of integration 


c,=— ( 1 - 3(5) 
1 2 

(B19) 

^2 = 3C 

(B20) 

-3C 4 


Cg (1 + cp ) 

(B21) 


where 

<p = l+ A = l+ - (B22) 

R 


The pressure distribution in the vapor gas is found by substituting 
equations (B17) and (B18) into the nondimensional form of equations (3) 
and (4) and solving for pressure. 
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j&fl.-jjljjfl.einc+ia i n (sec — \ + — 

C 4 C \ 2/ C 4 

q 

The constants C 4 and C can now be determined from the 
balance, equation (12), and the boundary condition, equation 

(p, 0) -* (1,0*) then p - p Q (B24) 

Applying condition (B24) to equation (B23) gives 


=^2. + 3(1 + <p) - 12 In (sec — \ (B25) 

C 4 C 4 V 2 ) 


(B23) 

static force 
( 11 ), 


The remaining constant C 4 may now be determined from equation (12), 
using equations (B7), (B8), (B15) to (B20), and (B25). 


W 


d 


C 4 = 


37TW* 

(cos^ 9* - 2 cos 0* + 1 )cp 


(B26) 


Therefore, utilizing Cp Cg, Cg, and C 4 , the velocity and pressure 
distributions are known functions of cp where <p - 1 + A. 

Next, the energy equation is considered. Afterwards, the solution to 
the energy and momentum equations are combined in the interface energy 
balance to obtain a solution for the heat-transfer coefficient, 

ENERGY EQUATION 
The solution of the energy equation (6) is 
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The temperature gradient at each surface becomes 


(B27) 


30 _ 

9? 'A 


(B28) 


30 1 

3 ? ?=<? 


(B29) 


INTERFACE ENERGY BALANCE 

The velocity and temperature distributions are now expressed in terms 
of two unknowns, A and 0*. The interface energy balance is used to deter- 
mine A in terms of 0*„ 

Substituting the values of constants Cj to C 3 from equations (B19) to 
(B21) into equation (B17) and evaluating v^ at cp gives the interface 
velocity 


JL = . 1/ 1 Z 3| ? + 3 - <p\ (B30) 

C 4 2 U 2 / 

where C^ is a function of 0*, cp, and the parameter W^/w* (see eq, (B26)). 

Nondimensionalizing equation (13) of the text, the nondimensional inter- 
face energy balance becomes 


v = -si® 
<P 


S=<P 


(B31) 
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where 


S = 


Pr 


C p (T d - T s> 


X* 


(B32) 


and X has been replaced by X* to accommodate some effects of convection. 

Substituting equations (B26), (B29), and (B30) into (B31) gives (for 

A « 1) 


a 4 /6rw*S\ / ^2 o „ 0 n* 

A = \ (cos 9* - 2 cos 9* + 1) 


W 


d 


' 


(B33) 


The problem is now reduced to finding the value of 6* from the bal- 
ance of forces on the interface, 

FORCES ON THE INTERFACE 

Within the supporting liquid, there are two forces to be considered, 
buoyancy and surface tension. Buoyancy is composed of two terms: 

(a) fluid displacement by the sphere to an angle of 0* and (b) fluid head 
as a result of fluid displacement by surface curvature. We now consider 
each of these forces separately. 

Surface tension support . - If a load is placed on an axisymmetric shell 
or membrane structure, it is supported at the circular interface or load 
circle. Here the effective load is supported by a component of surface 
tension 


F = 2uRo sin^ 9 * 
a 


(B34) 
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which may also be found by simply integrating the lift component of the sur- 
face tension pressure rise over the shell (i.e. , 




R.2 cos 9 sin 9 d0 d<& 


(B 3 5) 


Fluid displacement by the sphere. - The "pressure forces" acting on 
the spherical interface are illustrated in figure 3„ The buoyancy resulting 
from these forces is obtained by integration; however, one must note that 
0* and Z Q , while constants, are also unknown and must be determined 
from the force balance and interface curvature. 


B= J 0 



9 * 


0 


P^ cos 9 dA-^ 


'277 /~7T-0* 

'0 Jo 


P y cos 9 dAg 


(B36) 


Integrating equation (B37) gives the buoyancy in terms of 9* and Z Q as; 


F b = tt(p 


l 


p) Ji- R 2 <2R 
g n 


iP t +p) 


cos 0* , cos 0* 
+ — - 


3 (p t - p) 


+ Z 0 (l - cos'* 0*)} 


(B37) 


Fluid displacement by interface curvature. - As 9 increases beyond 
9*, the free surface is determined by solving the Young-Laplace capillary 
equation. 


D x u 


u 


X 


z 


(B38) 
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D v z = tan cp * 


u 

± Vl - u 2 


(B39) 


where 


and 


u = sin cp ^ 



L J 


(B40) 


where X is the horizontal coordinate. 

For 0* > 7t/2, tan cp > 0 and the + sign in equation (B39) is used; 
for 0* < w/2, tan cp < 0 and the sign in equation (B39) is negative. A 
fourth order Runga-Kutta routine was applied to equations (B38) and (B39) 
and values of 0* and z Q were determined for the case 


R = 0 . 3 cm 


P t - p = 0. 8 g/cc 
p g = 1.0 g/cc 
a = 8. 8 dyne/ cm 

which is shown as figure 7. While the forces are balanced to less than 

0. 3 percent, the exact solution was not obtained. For an exact solution, 

1 ° 

see reference 10. In this case 0* is approximately USt; (~1.98 rad) 
and the initial head is approximately 0. 148 centimeter. 
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A curve of 0* as a function of Bond number is given as figure 4, and 
can be used for approximating the value of 9 * for the most dense (heavi- 
est) levitated sphere. 

HEAT TRANSFER COEFFICIENT 
Drop Surface 

The Nusselt number for the sphere is given by 


h d R 

k 


3© 

a? 


C=i 


(B41) 


Substituting equations (B28) and (B33) into equation (B41) gives (assuming 
A « 1) 


h d R 


= 1 + 


W 


d 


a/4 


(B42) 


6ttw*S(cos z 0* - 2 cos 0* + 1) 


Liquid Interface 

The Nusselt number at the interface 


hjR 

k 


30 

a? 



(B43) 


Substituting equations (B24) and (B33) into equation (B43) gives 


hjR 

k 


1 + — = _l + 
A 


W 


d 


6ttw*S(cos^ I* - 2 cos |* + 1}_ 


1/4 


(B44) 
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Nondimensional Forms 

It can be shown that nondimensional forms identify the basic terms 
governing the solution. The parameter W d /w*S can be expressed as 


^d _ 4rr p d 
w*S 3 a p , - p 


(B45) 


Substituting equation (B45) into equations (B44) and (B42) indicates that heat 
transfer is influenced by the usual R and the interrelation between the 

ci 

density ratio [p^/ip^ ~ p)}> and Bond number. 


h^R 

k 

h^R 

k 


= 1 + 


= -1 + 


R k i 


9 a 


Pd \ 1 

Pi - PJ 

' p d \ r 

* - p ) f(e*) 


1/4 

1 1/4 


(B42a) 


(B44a) 


where 

f(0*) = cos 2 6* - 2 cos 9* + 1 


(B46) 
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TABLE IL - FREEZING DATA ON SMALL SPHERICAL WATER 
DROPS FLOATING ON LIQUID NITROGEN 
AT AMBIENT CONDITIONS 


Volume of 
drop, 
cm ^ 

Radius 

of 

drop, 

cm 

Time for 
dye 
change 
to be come 
visible, 
sec 

Time for 
drop to be 
completely 
frozen, 
sec 

Time for 
drop to 
fall to 
bottom , 
sec 

t Initial drop 
temperature , 
K 

a 0.09 


7 

16 

27 

298 

. 113 

a . 3 

8 

17 

23 

298 

a . 075 

.26 

9 

17 

28 

298 

a . 11 

.295 

9 

18 

27 

298 


a Measured parameter. 








r- LIQUID OR SOLID VOLUME 



Figure 1. - Schematic model for evaporation from a spherical interface. Figure 2. - Schematic of pressure and shear forces acting 

on the drop. 



Figured. - Submersion angle 0 on . for the optimum dens- 
ity ratio. 


Figure 3. - Forces acting at the vapor-liquid interface. 
Fluid displacement by the sphere. 




Figare 5. - Calculated surface temperature for a 0. 61 centimeter di- 
ameter water sphere on liquid nitrogen as a function of time. 


LEIDENFROST VALUES 

A SATURATED WATER ON STAINLESS STEEL 
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Figure 6. - Effect of surface on Leidenfrost temper- 
ature. 


APPROXIMATE INTERFACE SOLUTION 
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